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The conductance of disordered wires with symplectic symmetry is studied by a random- 
matrix approach. It has been beheved that Anderson locahzation inevitably arises in ordinary 
disordered wires. A counterexample is recently found in the systems with symplectic sym- 
metry, where one perfectly conducting channel is present even in the long-wire limit when 
the number of conducting channels is odd. This indicates that the odd-channel case is essen- 
tially different from the ordinary even-channel case. To study such differences, we derive the 
DMPK equation for transmission eigenvalues for both the even- and odd- channel cases. The 
behavior of dimensionless conductance is investigated on the basis of the resulting equation. 
In the short-wire regime, we find that the weak-antilocalization correction to the conduc- 
tance in the odd-channel case is equivalent to that in the even-channel case. We also find 
that the variance does not depend on whether the number of channels is even or odd. In the 
long-wire regime, it is shown that the dimensionless conductance in the even-channel case 
decays exponentially as (govcn) ~> with increasing system length, while (godd) ^ 1 in the 
odd-channel case. We evaluate the decay length for the even- and odd-channel cases and 
find a clear even-odd difference. These results indicate that the perfectly conducting channel 
induces clear even-odd differences in the long-wire regime. 
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1. Introduction 

The statistical aspect of phase-coherent electron transport in quasi-one-dimensional dis- 
ordered wires is characterized by three universality classes.^ The universahty classes describe 
transport properties which are independent of the microscopic details of disordered wires. The 
orthogonal class consists of systems having both time-reversal symmetry and spin-rotation in- 
variance, while the unitary class is characterized by the absence of time-reversal symmetry. 
The systems having time-reversal symmetry without spin-rotation invariance belong to the 
symplectic class. In the metallic regime where wire length L is much shorter than the local- 
ization length ^, the weak- localization (weak-antilocalization) effect arises in the orthogonal 
(symplectic) class, while the effect disappears in the unitary class. If the system length L 
satisfies -L 3> ^, the conductance decays exponentially with increasing L in all the classes. 
The mean-free path dependence of ^ is different from class to class. We can systematically 
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analyze such differences by using two different nonperturbative approaches. The first is the 
scaling approach of Dorokhov,^ and Mello, Pereyra and Kumar, ^ which is based on random- 
matrix theory.^ The second is the super- Fourier analysis by Zirnbauer,^'^ which is based on 
supersymmetric field (SSF) theory of Efetov.''' 

This work is concerned with the scaling approach. In terms of the Landauer formula, 
the dimensionless conductance g oi a disordered wire is given by the sum of transmission 
eigenvalues {T^} of individual channel a as 



where N is the number of conducting channels. The statistical property of the conductance is 
determined by the probability distribution function P{Ti,T2, ...,Ti\i; L) for the transmission 
eigenvalues. The scaling approach is based on an evolution equation of P(Ti, T2, Tjv; L) with 
increasing L, which is called the Dorokov-Mello-Pereyra-Kumar (DMPK) equation. This 
eqTiation was derived by Dorokhov for the unitary class, and by Mello, Pereyra and Kumar for 
the orthogonal class. Their derivation of the DMPK equation rests on the isotropy assumption 
of the distribution of scattering matrices. It is possible to replace the isotropy assumption by 
the weaker assTimption of equivalent scattering channels.^ Macedo and Chalker^ extended the 
DMPK equation to the symplectic class based on the isotropy assumption. 

Recent studies show that the transport property of the symplectic class is very unique 
compared with the other universality classes. ^'^"^^ From the scattering symmetry *r = — r 
generally holds in the symplectic class (r: reflection matrix), we can show that there exists one 
channel which is perfectly transmitting without backscattering when the number of conducting 
channels is odd. The presence of a perfectly conducting channel is first pointed out in the study 
of carbon nanotubes.^^"^^ Note that the perfectly conducting channel is present even in the 
limit of L ^ 00. This indicates the absence of Anderson localization in the odd-channel case. 
Since the odd-channel case is very different from the ordinary even-channel case, we must 
separately consider the two cases. That is, the symplectic class must be divided into the two 
subclasses: the symplectic class with an even number of channels and that with an odd number 
of channels. ^'^ Let Na {a = even or odd) be the number of conducting channels. We consider 
the dimensionless conductance gcvcn hi the even-channel case with Neven = 2m and (/odd in 
the odd-channel case with A'odd = 2m -|~ 1, where m is an integer. In the even-channel case, 
the number of independent transmission eigenvalues is m due to Kramers degeneracy, so the 
dimensionless conductance is expressed as gcvcn = '^'n^=i'^a- In the odd-channel case, one 
eigenvalue corresponding to the perfectly conducting channel is equal to unity, say T2m+i = Ij 
and other eigenvalues are doubly degenerate. Thus, we can express as (/odd = 1 + "^^1^=1 ^fl- 
it has been shown that 5even — > in the long-wire limit of L ^. In the odd-channel case, we 
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expect 

5odd ^1, (2) 

due to the presence of the perfectly conducting channel. Such a dramatic even-odd difference 
does not appear in the orthgonal and unitary classes. 

In this paper, we study the transport property of the symplectic class on the basis of the 
scaling approach. Our attention is focused on differences between the even- and odd-channel 
cases. We derive the DMPK equation for both the even- and odd-channel cases based on 
the assumption of equivalent scattering channels.^ In the even-channel case, we reproduce 
the result by Macedo and Chalker.^ Our assumption used in deriving the DMPK equation is 
weaker than the isotropy assumption of Macedo and Chalker. The DMPK equation for the 
odd-channel case is different from that for the even-channel case reflecting the presence of the 
perfectly conducting channel. We study the behavior of the dimensionless conductance on the 
basis of the resulting DMPK equation. The behavior qualitatively changes whether system 
length L is longer or shorter than conductance decay length (a = even or odd), where ^even 
is equivalent to the ordinary localization length. Thus, we separately consider the short- wire 
regime of L ^ and the long- wire regime of L ;» ^q, but treat the even- and odd-channel 
cases in a unified manner. 

In the short-wire regime, we obtain the ensemble average and variance of the dimension- 
less conductance using the perturbative treatment by Mello and Stone.^^ It is shown that the 
weak-antilocalization correction to ^odd is equivalent to that to g'even- Furthermore, we show 
that the difference between ^o^d and g'even is very small. This means that although one addi- 
tional channel is perfectly conducting in the odd-channel case, the even-odd difference does 
not clearly appear in the short- wire regime. The reason for this is explained from the view- 
point of the repulsion between transmission eigenvalues. We show that the repulsion from the 
perfectly conducting eigenvalue {T2m+i = 1) to other eigenvalues {Ti,T2, ■■■,Tjn) is essentially 
important. 

In the long-wire regime, we obtain the probability distribution function from the DMPK 
equation by using the treatment by Pichard.^^ On the basis of the resulting distribution 
function, we study the statistical property of the dimensionless conductance in both the even- 
and odd-channel cases. We show that the averaged ga behaves as (^even) — *• and (^odd) 1 
with increasing L. The behavior of ga in this regime is characterized by the conductance 
decay length ^q. To obtain ^q, we define Pq, as Peven = 5even/2 and Podd = (ffeven - l)/2, 
and evaluate the average of IuPq. The conductance decay length is obtained by identifying 
exp[(lnPQ,)] = exp[— 2L/^q;]- We find that ^odd for N^dd = 2m -|- 1 is much shorter than ^even 
for iVeven = 2m. This indicates that g^dd decays much faster than ^even with increasing L. 
Similar to the short-wire regime, the reason for this is explained from the viewpoint of the 
eigenvalue repulsion. Again, we show that the repulsion from T2m+i = 1 to other eigenvalues 
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is essentially important. 

Our results for both the short- and long-wire regimes are in qualitative agreement with 
those of the super- Fourier analysis.^' 

The present paper is a detailed and enlarged version of refs. 14 and 18. We here comment 
on major differences between the present paper and refs. 14 and 18. The derivation of the 
DMPK equation for the odd-channel case and the analysis of (^odd in the long-wire limit have 
been briefly described in ref. 14. However, a detailed comparison between the even- and odd- 
channel cases has not been presented there, so the even-odd difference originated from the 
perfectly conducting eigenvalue is less clarified than the present paper. In ref. 18, the average 
and variance of Qa in the short-wire regime have been evaluated in terms of a perturbation 
theory with the expansion parameter . The average {ga) is obtained up to 0(m^^) in ref. 
18, while the result up to 0(m^^) is given in the present paper. Thus, the present analysis 
for the short-wire regime goes one step beyond that in ref. 18. Several typographical errors in 
ref. 14 are corrected in the present paper. 

In the following, we restrict our consideration to the symplectic class. The outline of the 
present paper is as follows. In the next section, we first adapt the scattering approach to our 
problem, and then derive the DMPK equation for both the even- and odd-channel cases in a 
unified manner. The influence of the perfectly conducting eigenvalue becomes clear through 
its derivation. We also derive the scattering symmetry for the reflection matrix, from which 
the presence of the perfectly conducting channel is proved. In §3, we study the statistical 
properties of dinicnsionless conductance [a = even or odd) focusing on differences between 
the even- and odd-channel cases. We separately treat the short-wire regime of L <C and 
the long- wire regime of L ^ We obtain the average and variance of ga in the short- wire 
regime. A detail of the derivation is given in Appendix. In the long-wire regime, the average 
of In Tq, is obtained to evaluate the conductance decay length. We also obtain the average of 
ga and the variances of ga and InFo,. The influence of the perfectly conducting eigenvalue is 
discussed to explain the observed even-odd differences. In §4, we compare our results with 
those obtained from the SSF approach. The conclusion is presented in §5. 

2. Scattering Approach to Disordered Wires 

2.1 Scattering and transfer matrices 

To apply the scattering approach to our problem, we consider a disordered wire sandwiched 
between ideal leads as illustrated in Fig. 1. Let (a = even or odd) be the number of 
conducting channels. In this system, electron transport is described by the scattering matrix 
S, which is a 2Na x 2Na matrix given by 




(3) 
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l¥a->- 



l¥a-> 



Fig. 1. Disordered wire of length L sandwiched between two ideal leads. In the left lead, |V'o+) and 
I'^o-) correspond to the incoming and outgoing states, respectively. In the right lead, |V'a+) and 
\i}a-) correspond to the outgoing and incoming states, respectively. 



with Na X Na reflection matrices r and r' (reflection from left to left and from right to right) 
and transmission matrix t and t' (transmission from left to right and from right to left). 
The Hermitian matrix tt^ and i/i/^ have the common set of eigenvalues Ti,T2, ...,Tm in the 
even-channel case of A?even = 2m. Each eigenvalue is two-fold degenerate due to Kramers 
degeneracy. In the odd-channel case of N^dd = 2m -|- 1, there arises one additional eigenvalue 
which corresponds to the perfectly conducting channel. The dimensionless conductance is 
equal to the sum of all the transmission eigenvalues, so we obtain 

m 

5odd = l + 2^r„, (4) 

a=l 

m 

5even = 2 ^ ^ T^. (5) 

a=l 

We here derive the scattering symmetry for the reflection and transmission matrices in 
the symplectic class. We introduce the Hamiltonian H which is given by 

n = 'Ho + V, (6) 

where Ho describes the electronic states in the absence of impurities and V is the impurity 
potential in the disordered region. The above Hamiltonian has time-reversal symmetry 
without spin-rotation invar iance. This means that the ensemble of H is invariant under the 
transformation H UHU^^, where U is arbitrary symplectic matrix. In the absence of 
disorders (i.e., V = 0), we introduce the propagating wave |V'a+) the channel a (a = 
1,2,..., Na) from left to right (see Fig. 1). Note that |V'a+) is an eigenfunction of Hq. We 
introduce the time-reversal operator Ut given by 

Ut = KoC, (7) 

where C is the complex-conjugate operator and Kq is defined as 

Kq = -IGy. (8) 

In terms of the time-reversal operator, we define 

\'^a-)^UMa+) (9) 



5/28 



J. Phys. Soc. Jpn. Full Paper 

which represents the propagating wave in the channel a from right to left. Note that l^^-) 
is also an eigenfunction of Hq. Using eq. (9) and the relation olUf = —1, the eigenfunction 
Wa+) is given by 

\'4^a+) = -UMa-)- (10) 

Since our system belongs to the symplectic class, H. has the time reversal-symmetry as 

Utuu^^ = n. (11) 

To observe the symmetry of r, we define 

r,-„+ = (V'6-|r«|V'„+), (12) 

where is the T matrix given by 

T^ = V + V V + V V V-\ (13) 

E-Ho + ie E-Ho + ie E-H^ + ie ^ ' 

with E the energy of the scattering state. Prom eqs. (7), (8), (11) and (13), we find that 

UtT^U^^ = (T^)\ (14) 
Using eqs. (9), (10) and (14), the right-hand side of eq. (12) is rewritten as 

(Vb-|r«|V„+) = -(c/tV(,+ |T^|t^tV'„-) = -(V„-|r^lV'6+)- (15) 

Thus, from eqs. (12) and (15), we obtain rft-„+ = — r„-ft+. This immediately gives the relation 

na = -Tab- (16) 
To observe the symmetry of t and t', we define 

T,-,- = (^,-lT^|V;„-). (17) 
Prom eqs. (9), (10) and (14), the right-hand side of eq. (17) is rewritten as 

(V'b-lT^IV'a-) = (C/tV(,+ |r^|t^tV'a+) = (V'a+|r^lV';,+ )- (18) 

Thus, from eqs. (17) and (18), we obtain T^-a- = ^"0+6+- This immediately gives the relation 

t'ba = iab- (19) 

We observe that the transmission and reflection matrices satisfy-*^^ 

H = t', (20) 

V = -r and V = -r'. (21) 

Equation (20) holds for arbitrary systems with the time-reversal symmetry if an appropriate 
set of basisfunctions is chosen, while eq. (21) characterizes the peculiarity of the symplectic 
class. We can show with eq. (21) that there exists one channel which is perfectly transmitting 
without backscattering in the odd-channel case.^°'^^ 
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To study the evolution of transmission eigenvalues, it is convenient to introduce the trans- 
fer matrix M expressed by 



M = 

We employ the parameterization^' 

' e^ 



(1+7777^)^ 7? \ 

r?^ (l + r/'I'r/)^ / 

with 9 = ih, where h is an arbitrary Na x N^, Hermitian matrix and rj is an arbitrary x 
complex matrix. They satisfy 9^ = —9 and *77 = —rj. Equation (23) ensures both the flux 
conservation and the unique scattering symmetry eqs. (20) and (21). In terms of 9 and r), we 
can express the transmission and reflection matrices as 

t = e'^{l + r]rj^)~^, 

t' = (1 + ?7tr/)-ie-^*, 

r = —(1 + r]^r])~"^r]\ 

r' = e^r]{l + r]^r])~^e~^* . 

We introduce a unitary matrix v, reducing t't'^ to the diagonal form as 

vt't'K^ = diag{Ti,T2,...,TNj. 

Prom eq. (25), we derive t't'^ = (1 + T/^r?)"^ and find that 

vt't'K^ = (1 + r/1'r/)"\ 

where fj = v*r)v^ . Note that rj^rj has the eigenvalues Ai, A2, A;v„ which are given by 

1-Ta 



(22) 



M 







(23) 



(24) 
(25) 
(26) 
(27) 



A„ = 



Ta 



(28) 
(29) 
(30) 



Since fj is the A^^ x Na skew-symmetric matrix with *?} = —fj, we obtain 

'0 ^ \ 



V 



< 



(31) 



for the even-channel case with Neven = 2m, where ^ is an m x m diagonal matrix given by 

C = diag (VaT, v^, v^) . (32) 

In the odd-channel case with N^dd = 2m -|- 1, we assume without the loss of generality that 
the (2m -|- l)th channel corresponds to the perfectly conducting channel (i.e., X2m+i = 0). 
This means T2m+i = 1- We obtain 

/ e o\ 

??= ; • (33) 

V ... 0) 
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We see that 

■q^fj = diag (Ai, A2, A^, Ai, A2, A^) (34) 

for the even-channel case and 

fi^fj = diag (Ai, A2, Xm, Ai, A2, A^, 0) (35) 

for the odd-channel case. We shall use the notation, a = a + m iov 1 < a < m and a = a — m 
for m + 1 < a < 2m. Using a, we can rewrite the degeneracy relation = Ta+m (1 < a < "m) 
asTa = Ta{l<a< 2m). 
2.2 DMPK equation 

The distribution of the transmission eigenvalues depends on the wire length L. Here, we 
attach the small segment of 5x in the right-hand side of the system in Fig. 1 and discuss the 
resulting deviation of the transmission eigenvalues. Let t'-^ be the transmission matrix t' with 
the wire length L and t'2 be that with the wire length L + 6x. We express the transmission 
eigenvalues of the Hermitian matrices t[ti and tg^^ ^ {^a} and {Ta}, respectively. The 
transmission eigenvalues Ta and are related as 

fa = Ta + STa, (36) 

where 5Ta expresses a deviation. Wc focus on the eigenvalues from Ti to Tm and their distribu- 
tion which depends on L. It has been shown that the cvohition of the probability distribution 
function P(Ti, T2, T^', L) with increasing L is generally described by a Fokkcr-Planck equa- 
tion, which is usually called Dorokhov-Mello-Pcrcya-Kuniar (DMPK) equation in context of 
the scaling theory. -"^ In terms of {5Ta)ss and {5Ta6Ti,)ss, where (• • • )ss denotes the ensemble 
average with respect to the small segment of length 6x, the DMPK equation is given by^ 

To obtain an explicit form of the DMPK equation, we miist calculate the ensemble averages 
{5Ta)ss and {bTaST^ ss- To second order in peturbation theory, we obtain 

STa =Waa+ t -T ' 

where Wab is an element of the Hermitian matrix 

w = t'24 - i'A (39) 

in the basis where t'lt'l is diagonal. To relate t'l and t'2, we introduce the transfer matrix Mi for 
the wire of length L and the transfer matrix 5M for the small segment. The transfer matrix 
M2 for the combined system is given by 

M2 = 5MMi. (40) 
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We apply the parameterization given in eq. (23) to 5M. We shall take the weak-scattering 
limit, ^ where a moment higher than the second plays no role, so that SM can be expressed in 
the simple form 

/ 1 + r] \ 

SM=[ ^ • (41) 

Prom eqs. (22), (40) and (41), we can relate t'2 and t[ as 

t'2 = t[ (l-e*- + ^*^) . (42) 
Following Mello and Tomsovic,^ we assume that 

(^).. = {ifiss = 0, (43) 

{&abVcd)ss = {&abvld)ss = {VabVcd)ss = 0, (44) 

+ 7777^)^^ = 0, (45) 

{0abOcd)ss = Kl^,cd, (46) 

{eabe*Jss = l^lb,cd, (47) 

{Vabvld)ss = Klb,cd- (48) 

From eqs. (39) and (42), we decompose the Hermitian matrix w as 

U; = U;(1) +-U;(2)^ (49) 



where 



= vt[ (^6**2 + rj^r'y^f]^ t'lv^, (50) 
(2) = vt'i (^-r'^ri - 77VI) t'lvK (51) 



w 



We find that 



{W)ss = {W^'^)ss, (52) 

{ww),s = {w^^^w^^^),s. (53) 

We here take the weak-scattering limit ;^ Sx approaches and at the same time, we let scat- 
tering in the small segment become infinitely weak, so that 

<b,cd = lim (<5x)-i<fe ^rf (54) 

/t->0 

are finite quantities, while {Sx)~^ times a moment higher than the second vanishes. Accord- 
ingly, we calculate 

fa= lim{Sx)-^{STa)ss, (55) 

(5a;— »0 

K->0 

fab= lim{Sx)-^{STadn)ss, (56) 

ox—>0 

K-»0 
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instead of {5Ta)ss and {dTa5Th)ss- From eqs. (38), (52) and (53), we obtain 

fa = Um (fe)-i(^W).. + Hm (<5x)-i Yl ^^^("^f^ia (^7) 
lim(fe)-i(^(2>g)),,. (58) 

oa;— »u 

In accordance with 9^ = —9 and ^rj = —rj, we adopt the model^^ in which cr^^ is given by 
'^lb,cd = }^ini5xy^{9ab9cd)ss = -Sa,dSb,cCrab^ (59) 

^lh,cd = lim (5x)~^ {&abOcd)ss = Sa,cSb,d(^ab^ (60) 

f^a6,cd = Iim„(^a;)"^(r7a677*<i)ss = iSa,dSb,c - Sa,cSb,d)^ab- (61) 

ox— >0 

K->0 

Here, CTaf, represents the average reflection coefficient per unit length from b to a, while 
represents the average transmission coefficient. They satisfy X^^CTafe = Ylb^ab since (^^ + 
VV^)ss = 0. To proceed, we must assume that aab is expressed in a simple form. For the even- 
and odd-channel cases, we adopt the simplest choice^^ 

where / is the mean free path. This choice corresponds to the equivalent channel model, ^ and 
ensures the absence of backscattering (i.e., = 0)^°'^^ which characterizes the peculiarity 
of the symplectic class. Equation (62) leads to 

Y(^ab = Y(^'ab = J- (63) 



From eqs. (59)- (63), we obtain 



lim{Sx)-\w^^^)ss = -]{vt[t'^v%a + ^_^.^ («V)„, • 5^(1 - n) 



(64) 

^hm (fe)-l(^lJ^£)),, = ,^^_^^^ {{<Ar'A^'')aa ■ {v^A^'^)bb 

+ (Ki^^)aa • {vi!^rlr'Av'^)bb - {vArTA'v)ab 

X {y*i!^AAv'^)ba - {vi!^A%'v)ba ' {y*i!^Ai!lv\b} . (65) 

X {v%*r[t'lv^U - (vtWH'vha ■ {v%*r[t'lv^)ba}. (66) 



10/28 



J. Phys. Soc. Jpn. 

Equations (24)-(27) give 

vt\t'lv^ = {l + fi\fii)-\ 
vt[r'^r[t'^v^ = (1 + 77|^i)~^f7|r7i, 

vt[r'^%^v = {1 + rilmym 
vXr[t'lv^ = mil + rjlfii)-^, 
where rji = v*r]ivK Prom eqs. (34), (35) and (57)-(70), we obtain 
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(67) 
(68) 
(69) 
(70) 



1 T 



= -TTa + 



\ 



V 



6=1_ 



Tg + Tb- 2Tan 

Ta-n 



and 



/ " ' (TV, - 1)1 



fab 



'^-Ta- 5a,odd + ^ ^ 



V 



1 



b=l_ 
6(^a,a) 



j{Sa,b + S-a,b)2T^{l-Ta 



(Na - 1)1 

for 1 < a, 6 < 2m. For a = 2m + 1 in the odd-channel case, we also obtain 



(71) 



(72) 



/2m+l = 0, (73) 
/2m+l b = 0. (74) 

In cq. (71), we have introduced ^Q^odd which is defined as (Jo^odd = 1 for a = odd and for 
a = even. It should be emphasized that S^^odd indicates the contribution from the perfectly 
conducting eigenvalue. Substituting eqs. (71) and (72) into eq. (37), we obtain 

a=l V ^ / 

where s = L/l is the normalized system length. Changing variables from to Aa = (1 — 
Ta)/Ta^ we finally obtain the DMPK equation 



dP{\i,...,\m\s] 

ds 



-^^i^^dxA^^^^^'^^'^'^dK 

a=l ^ 



d f P{Xi,...,Xjn;s) 



(76) 



The even-odd difference is described by J,, 

m—l m 

Jeven = J] 11 " l"' C^^) 

6=1 a=b+l 

m m—l m 

Jodd=n-^cx n n (^s) 

c=l 6=1 a=6+l 

Note that the factor H^i '^c <^dd represents the eigenvalue repulsion to Ai,A2,-- - , A^ 
arising from the perfectly conducting eigenvalue of A2rn+i = 0- This eigenvalue repulsion is 
the origin of all the even-odd difference. 
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3. Transport Properties 

To study the transport properties in disordered wires, it is convenient to introduce Fq, 
{a = even or odd) which is defined as 

m m -. 

a=l a=l ^ ^ 

Prom eqs. (4), (5) and (79), the dimensionless conductance for the even- and odd-channel 
cases is given by 

9a = Sa,odd + 2ra. (80) 

3.1 Short-wire regime 

We adopt MeUo and Stone's method^^ to obtain the average and variance of ga, 

{9a)=SaMd + '^(^a), (81) 

Varb«]=4((r2)-(r„)2), (82) 

where 

„ in 

(•••)= / l\d\a---P{Xl,...,Xm; S). (83) 

a=l 

Prom eq. (76), we can derive the evolution equation for the ensemble average {F{Xi, Xm)) of 
an arbitrary function F(Ai, Am).^^ Multiplying both sides of eq. (76) by F and integrating 
over {Xa}, we obtain. 



a=l 6=1 ^« - 



+<5a,odd^f;{2(l + A„)|£}y (84) 



Let us derive the evolution equation for (Pa). Setting F = Pq, in eq. (84), we obtain 

(iV« - 1)^1^ = (-2P2 + P,,2 - 2Pa,5«,odd>, (85) 

where V^q was defined as P^g = Fa - To solve eq. (85), we need the quantities (P^) and 

(Pa2)- So we set F = (Pa2) in eq. (84) and obtain 

(iV„ - 1) = (4P2 - 8P„Pa2 - 2(1 + 25c«,odd)ra2 + 4Pa3> . (86) 

New quantities (P^), (PaPa2) and (Pas) appear in the right-hand side of eq. (86). We will 
observe that every time we write down an evolution equation, we find, on the right-hand 
side, new quantities that have not appeared before. This means that we cannot obtain a 
closed set of coupled equations, and thereby cannot obtain an exact solution. To overcome 
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this difficulty, we employ the method by Mcllo and Stone^"*^ to obtain an approximate solution 
for (Fq,) and (F^). This method is based on a series expansion with respect to m^^, and is 
applicable to the short-wire regime oi s <^ m. Thus, we assume in our following argument that 
m is much larger than unity. We derive the evolution equation for (Ta), as well as those for 
the ensemble averages {T^Ta2), i^a'^as), (rSr^2)> i^a'^ai), a2^a3) and (r^r^g)- Setting 
F = rS, T^Ta2, rgr^s, r^r^g, V^TaA, ^^^a2^a3 and rSr^2 in eq. (84), we obtain the 
evolution equations 



(Na-l 



on) 

ds 



{Na - 1) 



ds 



{Na - 1) 



a(rSr^4) 

ds 



^ " ^ ds 



{-2pTl+'+pTP,-'Ta2 

+p{p - i)rr'(ra2 - r«3) - 2pr^(5,,odd>, (87) 

{4Tl+'-2ip + A)TP+'ra2 

-2{i + {p + 2)(5a,odd}r^r«2 + 4r^r«3 
+prp,-'Tl^ + 4prr'(r«3 - r«4) 
+p{p-i)rp-^{rl,-ra2ra3)), (88) 
{-2{p + 6)r^ir„3 + i2rrir„2 - er^^r^^ 
-2{3 + (p + 3)(5„,odd}rsr,3 + 9rpr„4 
+prr^ra2ra3 + 6prr\r«4 - r^s) 
+p(p-i)rr'(r„2r„3-r23)), (89) 

(-2(p + 8)rs+ir22 + 8rs+2r,,2 + ^r^^ 
-2{2 + {p + 4)5e«,odd}rpr22 + 8rsr,,2ra3 

- i)rr'(ra2 - r„3)r22 + 8rp (r«4 - r^s) 
+8prr^r<,2(ra3-ra4)>, (90) 
{p{p - i)rp-'(ra2 - r«3)rc4 + prr'r„2r«4 
~2{p + 8)rp+ir„4 - 2{(p + 4)<^,,odd + 6}rprc,4 
+i6rf;r«5 + 8prp^-\ra5 - r^e) + i6rg+^r,,3 
+8rsr22-i6rpr«2r«3>, (9i) 

{p{p - i)rf;-2rc,2r«3(r«2 - rc.3) 
+4prr'r„3(ra3 - ra4) + 6prr'rc,2(ra4 - r^s) 
+i2rg(r„5 - r„6) + 4rpr2 3 + prr'r2 2r„3 
-2{4 + ip + 5)(5„,odd}rsrc,2ra3 + 9rpr<,2ra4 
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-er^r^^), (92) 

+24rsr«2(ra4 - r^s) - 2p{p + i2)tp,-^'tI^ 
+i2rpr22r„3 + i2prr'r^2(r«3 - r^^) 

+l2TP+^Tl^ - 2{3 + (p + 6)5„,odd}rSr3 2>. (93) 

We seek the solution for the above equations as a series in decreasing powers of m. To do so, 
we expand the averages as 

(r^;) = mPU,p,oi^) + mP-'U,p,ii^) + mP-^U,p,2is) + mP-^U,p,3i^) + ■■■, (94) 

(rPr«2) = mP+'ja,p+iAs)+mPja,p+i,iis) + mP-'ja,p+i,2{s) + --- , (95) 

(r^r^a) = mP+^ka,p+i,o{s) + mPka,p+i,i{s) + --- , (96) 

(r^r^^) = mP+X,p+2,o{s)+mP+X,p+2,i{') + --- , (97) 

(rSr„4) = mP+Ha,p+i,o{s) + -- - , (98) 

(rsrc,2ra3) = mf+V,p+2,o(s) + • ■ ■ , (99) 
(r^r^2) = m^'+V,p+3,o(s) + ---, (100) 

where fa,p (s) and ya,p,n{s) are functions of s 

and satisfy the initial conditions 

fa,p,n{0) — Ja,p,n(0) = ^a,p,n(0) = ^a,p,n(0) = ia,p,n(0) = ^^a,p,n(0) = J/a,p,n(0) = (^n,0- (101) 

Prom eqs. (87)-(101), we evaluate fa,p,o{s), fa,p,i{s) ■ ■ ■ ■ The detail of calculations is given in 
Appendix. Substituting eqs. (A-8), (A-15), (A-19) and (A-20) into eq. (94), we finally obtain 



+ rnP-^ P 



360(1 + s)P+4 



+ m 



{lip - 9 + (15p - 15)Sa,odd}s'' 

+ {36p - 24 + (120p - 120)6^,odd}s^ 
+ {90p - 60 + (240p - 240)(5a,odd}s^ 
+ {120p - 150 + (180p - 180)(5a,odd}s^ 
+ {90p - 90 + (45p - 45)(5c,,odd}s'^ 
{(16V - 369p + 196) 



45360(1 + s)P+6 

- (693p2 - 1953p + 1260)(5a,odd}s^ 
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+ {(756p2 - ISmp + 630) - (554V - 15624p + W206)6a,odd}s^ 

+ {(1890p2 - 2448p + 1008) - (19404^^ - 54684p + 36288)(5a,odd}s'^ 

+ {(2520p^ - 2142p + 1260) - (36288p^ - 106218p + 73458)(5a,odd}s^ 

+ {(1890j9^ + 3402^) - 4914) - (39690j92 - 120960p + 83916)5a,odd}s^ 

+ {(9450p - 24570) - (24570^^ - 75600p + 49140)(5a,odd}s'' 

+ {(18900p - 24570) - (6615^^ - 19845p + 11340)(5a,odd}s^ 

+ (5670J9 - 5670)5^ + • ■ ■ . (102) 

Setting p = 1 in eq. (102) and substituting the resulting expression for (r^) into eq. (81), we 
obtain 

- -=(25^ + 18/ + 93s^ + 315^^ + 378s^ + 2205s^ + 630s^) 

3780(1 + s)^^ ^ 

+ O (m-3) , (103) 

1 s3 ,1 

1 1- m 

1 + s 3(1 + s)3 90(1 + s)5 



(5even) = 2m:p^ + + m- \^^^ , ^,, (/ + 6s^ + 15s^ - 15s=^) 

1 



- -=(25^ - 3/ - 75s^ - 273s^ - 63s^ + 2520s^ + 9455^) 

3780(1 + s)7^ ^ 

+ O (m-3) . (104) 

Now, we discuss the behavior of the averaged conductance in the short-wire regime of 
s <^m. The first term in (ga) represents the classical contribution. The second term in {ga) 
represents the weak-antilocalization correction. We observe that the weak-antilocahzation 
correction to (^odd) is equivalent to that to (geven)- We also observe that the third term in 
(ffodd) is equivalent to that in (^even)- We find a small even-odd difference in the forth terms of 
(ffodd) a'Hd (geven), although their leading order corrections are equivalent to each other. This 
indicates that, in the short-wire regime, a notable even-odd difference does not appear in the 
averaged conductance. We here discuss the reason why (^odd) is almost equivalent to (^even)- 
Note that the perfectly conducting channel is present only in the odd-channel case, and that 
if we neglect this special channel, the total number of conducting channels in the odd-channel 
case is equivalent to that in the even-channel case. Thus, one may expect 

(5odd) = (5even) + 1- (105) 

However, contrary to this hypothesis, our result indicates that ((7odd) ~ {dcvcn)- To explain 
this, wc must consider the repulsion acting between transmission eigenvalues of {Ta}. In the 
odd-channel case, one transmission eigenvalue T2m+i is equal to unity due to the presence 
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of the perfectly conducting channel. Although it positively contributes to (godd), the other 
eigenvalues {Ti, T„j} are reduced due to the repulsion from it. To observe this, we set p = 1 
in eq. (102) and present the expression for Tq, 

(ra) = ^ (TT^ + ^(Y^ - ^(T^^"-odd + o (m-^) . (i06) 

The third term in eq. (106) represents this reduction. Our result indicates that the contribution 
to (^odd) from the perfectly conducting channel is almost canceled out by the reduction of the 
contribution from the other channels. 

Next, we obtain the variance Varf^a] for the even- and odd-channel cases. Prom eq. (102), 
we can evaluate (F^) and (Fq)^. Substituting these quantities into eq. (82), we obtain 

Var [godd] = Var [^even] 



15(1 + s)6 



(s*^ + + 15s* + 20s'^ + 15s^) 



315(1 + s)^ 

+ 504s'' + 525s'^ + 1050s^ + 315s2) + O (m'^) . (107) 

The variance does not depend on whether the total number of channels is even or odd. 
3.2 Long-wire regime 

We adapt Pichard's treatment^^ to our problem to study the behavior of ga in the long- 
wire regime. It is convenient to introduce a new set of variables Xa, related to Xa by^ 

Xa = smh^Xa, (108) 

where Xa > 0. In the odd-channel case, the perfectly conducting channel is characterized 
by X2m+i = 0. If we make a change of variables from Xa to Xa, the evolution equation for 
P{xi, ...,Xm', s) becomes 

dP _ 1 ^ d fdP d9^\ 

'd^ ~ 4(iV„ - 1) A. l^a^ + -Q^J ( ) 

with 

fia = - In ^ sinh 2xa ■ ■ (110) 
The explicit form of fi^ is given by 

m m—1 m 

= — ^ In I sinh 2xa\ — 4 ^ ^ In | sinh^ Xa — sinh^ Xb| 

a=l 6=1 a=6+l 

m 

- 2^1n|sinh^Xa|(5a,odd- (m) 

a=l 

In the limit of s/N^ = L/Nal — oo, we expect that the variables Xa{^ < a < m) are much 
larger than unity and are widely separated. We thus assume that 1 <^ xi <^ X2 ■ ■ ■ "^i Xm- 
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Under this assumption, we can approximate that $7^ ~ —2 ^^j^(4a — 3+2(5cK,odd)iCa+constant. 
Substituting this into eq. (109), we obtain 

where 

7„ = 2(Ar„ - 1). (113) 

The solution of the above equation is 

m rn ^ , I \^ 

P^{x,,...,xm;s) = (^) ' He-^y'^'-^V , (114) 

where 



V27rs/ 

a=l 



^"'""4a-3+W,odd- ^^^^^ 
We observe from eq. (114) that each eigenvalue obeys a Gaussian and that the even-odd 
difference appears in the distribution function. Prom the resulting distribution, we find that 
the root-variance -y/Var[x^ ~ \/sJla is indeed much smaller than their spacing 

{Xa+l - Xa) ~ s/7q- (116) 

Note that [xi) ~ s/7a. Thus, when 5/7^ ^ 1, the dimensionless conductance Qa = Sa,odd+'^^a 
is dominated by xi. That is to say, considering eq. (79) and relation of Ta = 1/ cosh^Xa, 
behaves as 

ra!«4e-2^i. (117) 
Equation (114) gives the distribution function for xi as 

We will discuss this cvcn-odd difference later. Using the above distribution function, the 
ensemble averages (xi) and (xf) are obtained as 

(m) = — (l + 2,5,,odd), (119) 

7a 

{xj) = A + + (120) 

7a 7a 

Based on eq. (117), we evaluate the ensemble-average — (Inra) = 2(xi). Prom eqs. (113) and 
(119), we obtain 

3L 
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We estimate the decay length ^q, of the conductance by identifying 

exp[(lnr«)] = exp[-2L/^„]- (123) 

We find that 

Codd = ^(^odd - 1)/ = ^ml, (124) 

Ceven = 2(Are^en - 1)/ = 2(2m - 1)1. (125) 

Note that ^odd < '^cven for any positive integer m. This indicates that ^odd decays much faster 
than geven- Prom eqs. (117) and (120)-(122), we also obtain 

VVax[lnrodd] = (126) 

VVar[lnre.en] = jj^^. (127) 



Note that y^VarpnTo,] << Klnr^)! for the even- and odd-channel cases. Next, we obtain the 
ensemble average (r^) as 

7q, / {1 + 2(5„,odd)s 
Tri^i — 



(Ta) cx dxi4e-2^i ■ exp 



la 

oce"«^, (128) 
where ^^dd ^nd ^'^^^^ are given by 

Cdd = (^odd - 1)^ = 2mZ, (129) 

Cen = 8(Are,e„ - 1)/ = 8(2m - (130) 

The relation ^'^^^ < ^'^^^^ holds in any positive integer m. The above results indicate that 

(5even) ~e-2VCe„, (131) 

(-^^odd) = (5odd) - 1 ~ e-2^/Cd . (132) 

We also obtain the variance of Fq, as VarfFeven] ~ e~"^^/^«^en and Var[rodd] ~ e~^^/^^odd. This 
results in 

Varbeven] ~e-2^/€-n, (133) 

Var[5odd]-^e-9V<dd. (134) 

It should be noted that in the limit of L/^^ ^ oo, (Fodd) and (Feven) are much smaller than 
V^Var[Fodd] ~ e^^-^/^^odd and y^Var [Feven] ~ e~^/^'<'^<"^ , respectively. 

As noted below eq. (76), the even-odd differences obtained above should be attributed to 
the eigenvalue repulsion from the perfectly conducting eigenvalue. The distribution of trans- 
mission eigenvalues has been given in eq. (114). We find that the distribution of {xi,X2, Xm} 
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in the odd-channel case shift toward larger values than those in the even-channel case. In- 
deed, we find from eq. (119) that the peak of Poddixi, 1) is at (xi) = 3s/7odd) while that of 
Pevenixi, 1) is at (xi) = s/7cvcin where 7odd ~ 7cvcn for large m. Since Xa and are related by 
Ta = 1/ cosh^Xa, the transmission eigenvalues {Ti,T2, ...,Tm} are reduced in the odd-channel 
case compared with the even-channel case. This results in the suppression of Todd- This ten- 
dency should be interpreted as a manifestation of the repulsion from the perfectly conducting 
eigenvalue. 

4. Comparison with the SSF Approach 

In this section, we compare our results for {ga) and Var[^cK] (a = even or odd) given 
in the previous section with those obtained from the supersymmetric field (SSF) approach. 
We first present the analytic expressions for (ga) and {g'^) using the Zirnbauer's super-Fourier 
analysis.^' ^'^^ These expressions have been derived previously except for {g^^^), so we describe 
their derivation very briefly. The super-Fourier analysis is based on the nonlinear a model 
derived by the SSF approach^^ and is valid for the thick- wire limit which is deflned by AT ^ oo, 
L/l oc with a fixed Nl/L, where N and I are the number of conducting channels and 
the mean free path, respectively. The conductance and its moments are expressed by the 
propagator, called the heat kernel, for the nonlinear a model. In the super-Fourier analysis, 
we expand the heat kernel in terms of the eigenfunctions of the Laplacian defined in the a- 
model space. Note that in addition to normal modes, the eigenfunctions contain a zero mode 
and its subsidiary series.^ Using the resulting expansion, Zirnbauer obtained the expression 
of the average conductance, 

(g) = iT(i, 1,1)+ I {^(i' n,n-2)+ T{i, n - 2, n)} 



2 

n=3,5,7,--- 



+ 2^ E 

Tii,n2=l,3,5,- 



POO / 

J dAA(A^ -I- 1) tanh ( J "'i"'2Pi (A, ni , 712) 



X Yl (-l + icrA-Fcrmi -Fc72n2)~^r(A,ni,n2) (135) 

0-,(Tl,(T2=±l 



with 



T{X,ni,n2) = exp | - (A^ + n? + ni - 1) ^1 , 



(136) 



pi{X,ni,n2) = + nl + nl-l. (137) 

Here, ^ in eq.(136) is a typical length scale characterizing the decay of the conductance. 
The first (second) term of eq. (135) corresponds to the contribution from the zero mode 
(the subsidiary modes), and the third term comes from the normal modes. Equation (135) 
indicates that (g) — 1/2 when L 00.^'^ This apparently contradicts the scaling theory, 
which gives (geven) — and (^odd) ^ 1- Note that the even-odd difference is not reflected in 
the Zirnbauer's argument. Brouwer and Frahm^^ introduced a parity operation and classifled 
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the eigenfunctions of the Laplacian into even-parity modes and odd-parity modes. The zero 
mode and the subsidiary modes have odd parity, so that the first and second terms are classified 
into the contribution from the odd-parity modes. The part of the third term with ni+n2 = 2 
(mod 4) is the contribution from the even-parity modes, while that with ni+n2 = (mod 4) 
is the contribution from the odd-parity modes. They pointed out that the odd-parity modes 
are unphysical in the even-channel case. We must neglect the contribution from the odd-parity 
modes in calculating (gcvcn), and double the contribution from the even-parity modes to take 
Kramers degeneracy into account. We obtain^^ 

(ffeven) = 2^ ^ / dAA(A^ -|- 1) tanh [ ^ j nin2Pi{X, ni, 722) 

ni,n2=l,3,5,- V / 

ni+n2S2 (mod 4) 

X Yl (-l + i(^A-Fo-iTii + cr2n2)"^r(A,ni,n2). (138) 

cr,<Ti,(T2=±l 

Although Brouwer and Prahm correctly pointed out that the odd-parity contribution must 
be excluded in the even-channel case, they were not aware that the odd-parity modes become 
essential in the odd-channel case. One of the present authors^^ pointed out that we must 
exclude the even-parity modes in the odd-channel case. Neglecting the contribution from the 
even-parity modes and doubling the contribution from the odd-parity modes to take Kramers 
degeneracy into account, we obtain^^ 

(ffodd) = r(i, 1, 1) + {Tii,n,n-2) + T{i,n-2,n)} 

n=3,5,7,--- 



+ 2^ V r dAA(A2 + 1) tanh (^] nin2Pi (A, 



ni,n2 



ni,n2=l,3,5, 
"■1+^12=0 (mod 4) 



X Yl (-1 + i(^A-Fo-ini -F(T2n2)~^r(A,ni,n2). (139) 

We observe that {g) = ((5cvcn) + (fi'odd))/2. On the basis of eqs. (138) and (139), we numerically 
calculate (geven) and (godd)- The results are displayed in Figs. 2 and 3. We observe that 
(5odd) = (5even) wheu L/^ <C 1 and that (^odd) ^ 1 and (gcvcn) ^ when L/^ » 1. 
The behavior of is in agreement with our scaling theory. Furthermore, we observe from 
eqs. (138) and (139) that (gcvcn) ~ e'^/^^ and (5godd) = (godd) - 1 ~ e~^^/«. That is, the 
characteristic length scale for (^godd) is a factor 8 shorter than that for (gcvcn)- This result is 
also in agreement with the scaling theory which yields Ccvcn/'^odd ~^ ^ large ni. Let us next 
consider the variance defined by Var[ga] = (ga) ~ ida)'^- We need the second moment (g^) to 
obtain Var[gQ]. Zirnbauer and co-workers^ found the expression for the second moment, 

(5') = ^r(i,l,l)+ Y Ul{n-lf + l){T{i,n,n-2)+T{i,n-2,n)} 
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Fig. 2. The averaged dimensionless conductance in the odd-channel case, {godd), and that in the 
even-channel case, (geven), as a function of L/^. 




Fig. 3. (fl'odd), (ffeven) and {6godd) = (ffodd) - 1 as a function of L/^. 



with 



+ 2* V / dAA(A2 + l)tanh ) nin2P2(A,ni,ra2) 

ni,n2=l,3,5,--^0 

X (-1 + io-A + (Jini + (72n2)"^T(A,ni,ra2) 

<T,<Tl,iT2=±l 



P2(A,ni,n2) = ^{nt + ra^ + (3A2 + l)(n? + ni) + 2A^-2A2-2} 



(140) 



(141) 



Again, the correct result for the even-channel case is obtained if we neglect the odd-parity 
contribution and double the even-parity contribution. The result is^^ 

(5even) = 2^ Yl / dAA(A2 + 1) tanh ( — j nin2P2(A, ni, 712) 

ni,n2=l,3,5,- ^ ^ 

ni+n2=2 (mod 4) 
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X Yl (-1 + it^A + aini + c72n2)"^T(A,ni,n2). (142) 

cr,CTi,(T2=±l 

Similar to the derivation of (^even)' ^® obtain 

(5odd) = r(i, 1,1)+ - + 1)) {^(i' n,n-2) + T(i, n - 2, n)} 

^ POO /^^\ 

+ 2^ / dAA(A2 + l)tanh ( — j nin2P2(A,ni,n2) 

ni,n2=l,3,5,--- ^ ^ 

»ii+»i2=0 (mod 4) 

X Yl (-l + i(^A + c7ini + c72n2)~^r(A,ni,n2). (143) 

0-,CTl,(T2=±l 

On the basis of eqs. (138), (139), (142) and (143), we numerically calculate Varf^even] and 




Fig. 4. The variance of the dimensionlcss conductance in the odd-channel case, Var[(7odd], and that 
in the even-channel case, Var[5even], as a function of L/^. 



Yar[go(\d]- The results are displayed in Figs. 4 and 5. We observe that Var[(7odd] = Var[(7cven] 
when L is much shorter than ^. This is in agreement with our analysis for the short-wire 
regime. With increasing L/^, Var[5(odd] very rapidly decreases to zero, while the decrease of 
Var[5(cvcn] is much slow. Furthermore, we observe from eqs. (138), (139), (142) and (143) that 
Var[5(cvcn] ^ e~^/^^ and Var[(7odd] ~ e^^^^'^^. Again, these results are in agreement with the 
scaling theory. 

5. Conclusion 

We have studied the quantum electron transport in disordered wires with symplectic 
symmetry. In the symplectic class, the reflection matrix r has the scattering symmetry *r = — r, 
from which the presence of a perfectly conducting channel without backscattering can be 
proved for the odd-channel case. The perfectly conducting channel is present even in the 
long-wire limit. This indicates the absence of Anderson localization in the odd-channel case. 
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Fig. 5. Var[5odd] and Varfgieven] as a function of L/^. 



Since the odd-channel case is very different from the ordinary even-channel case, we separately 
consider the two cases. To study the two cases within a unified manner, we adopt the scaling 
approach based on a random-matrix theory. The transport property is fully described by the 
probability distribution function of transmission eigenvalues. For both the even- and odd- 
channel cases, we have derived the DMPK equation which describes the evolution of the 
distribution function with increasing system length. In the even-channel case, the resulting 
DMPK equation is equivalent to the previous result for the symplectic class. The DMPK 
equation for the odd-channel case is different from that for the even-channel case reflecting 
the presence of a perfectly conducting channel. Using the resulting DMPK equations for the 
even- and odd-channel cases, we study the behavior of the dimensionless conductance for the 
short- wire regime of L/^q, <S 1 and the long- wire regime of L/^^ Ij where L is wire length 
and is the conductance decay length {a = even or odd). 

In the short-wire regime, we have derived the evolution equations for various quantities, 
such as (Ta) and {TaTci2), from which we obtain (ga) and Var[giQ,]. We have solved the resulting 
equations by using the series expansion with respect to s/m?^ We observed that the weak- 
antilocalization correction to (^odd) is equivalent to that to (^even)- Furthermore, we obtained 
higher order terms of the series expansion with respect to s/m and found a small even-odd 
difference between (^odd) ^^nd (^even)- This means that (^odd) ~ (^even)- We have discussed 
this behavior from the viewpoint of the eigenvalue repulsion from the perfectly conducting 
channel. We have also shown that Var[5odd] = Varf^even]- We thus concluded that clear even- 
odd differences do not appear in the transport properties in the short-wire regime. In the long- 
wire regime, we approximately solved the DMPK equation and studied transport property 
for the even- and odd-channel cases. We have shown that the dimensionless conductance in 
the even-channel case decays exponentially as (^even) — with increasing system length, 
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while ((7odd) ^ 1 in the odd-channel case. We evaluated the conductance decay length S,a by 
identifying exp[(lnrQ,)] = exp[-2L/^a], where = {g^ - 5a,odd)/2, and found ^odd < Ceven- 
That is, (7odd ~ 1 decays much faster than (/even- We have discussed this behavior from the 



viewpoint of the eigenvalue repulsion. We also obtained variance -^Var[ln r^] and observed a 
clear even-odd difference. 

From what has been discussed above, we conclude that the perfectly conducting channel 
induces the clear even-odd difference in the long-wire regime. 
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Appendix: Solution of the coupled equations 

We substitute eqs. (94)-(100) into eqs. (87)-(93) and equate the coefficients of various 
powers of m. Accordingly, we obtain the following differential equations, 

/a,p,o(s)+P/a,P+l,o(s) =0, (A-1) 
Ja,p,o(s) + iP + 3)ja,p+l,o(s) = 2fa,p+l,ois), (A-2) 
/4,p,l(s) + P/a,p+l,l(s) = ^ [Pja,p,0{s) + fL,p,o{s)da,even - 2p/a,p,o(s)5a,odd] , (A-3) 
4,p,o(s) + iP + 6)/a,p+l,o(s) = 4jc,,p+i,o(s), (A-4) 
K,p,oi^) + iP + 5)/Ca,p+l,o(s) = 6ja,p+l,ois) - 3la,p+l,ois) , (A-5) 
Ja,p,l(s) + iP + 3)ja,p+l,l(s) = 2fa,p+l,lis) + 2ka,p,ois) + ^(p - l)Za,o(s) 

+ \j'a,p,oi^)Keven - {I + {p + l)(^a,odd}ja,p,o(s), (A-6) 
/a,p,2(s) +P/a,p+l,2(s) = ^pja,p,lis) + ^p{p - l)ja,p-l,ois) 

- \p{P - l)fca,p-l,o(s) + ^/a,p,l(s)(^a,even 

-pfa,p,l{s)Sa,oM- (A-7) 
Note that the above equations arc closed for fa,p,o, ja,p,o, fa,p,i, ^a,p,o, ^a,p,0) 3a,p,i and fa,p,2- 
To obtain fa,p,3{s), we need equations for ya,pfi, Ua,pfi, ta,pfl, la,p,i, ka,p,i, ja,p,2 and fa,p,3{s) 
in addition to eqs. (A-l)-(A-7). By considering the series expansion of fa,p,o{s) at s = on 
the basis of eq. (A-1), we expect that 

We can easily show that eq. (A-8) satisfies eq. (A-1) for arbitrary s. Since the resulting power 
series for the other coefficients ja,p,o{s), fa,p,i{s), ■ ■ ■ at s = are not easily identifiable, it is 
difficult to derive all the coefficients in the above manner. 
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From the argument by Mello and Stone, a solution of coupled differential equations, 
eqs. (A-l)-(A-7), is given by the combination of polynomials in s divided by a power of (1 + s). 
Such a polynomial is written as 

= (i + g)p+n (^^^'' + (^r-is""^ + ■■■ + ao), (A-9) 

where n and r are integers and {a^} are coefficients which are independent of p. Regardless 
of integers n and r, Ap{s) satisfies the relation of 

Noting this observation, we can systematically obtain the solutions of eqs. (A-l)-(A-7). Indeed, 
if fa,p,o{s) satisfies eq. (A-10), we can rewrite eq. (A-1) as 

f'a,p,ois) + Y^Ja,pfi{^) = (A-ll) 

with the initial condition /a,p,o(0) = 1- Its solution is equivalent to eq. (A-8). 

We assume that ja,p,ois) satisfies eq. (A-10), and rewrite eq. (A-2) as a first-order linear 
differential equation. The solution is 

ja,p,o{^) = ^pT^^(2^' + + 6« + 3)- (A-12) 
Substituting eq. (A-8), its derivative and eq.(A-12) into eq. (A-3), we obtain 

/a,p,l(«) + Pfa,p+l,l{s) = ^ g^p+3 [2(1 - 3(5a,odd)s^ + 3(1 - 5(5a,odd)s^ 

-12s<5a,odd - 35Q,,odd] - (A- 13) 

We assume that fa,p,i{s) is expressed as 

/a,p,i(s) = P ■ Xa,p{s), (A-14) 

where Xa,p{s) satisfies eq. (A-10). Under the above assumption, we substitute fa^p^i{s) and 
fa,p+i,i{s) into the left-hand side of eq. (A-3) and obtain a first-order linear differential equa- 
tion for Xa,p{s)- We solve the resulting equation and obtain 

fa,p,l{s) = _|.^g)p+2 ((1 ~ 3(^a,odd)s^ - 6s^(5o,,odd - 3s(^a,odd) - (AT5) 
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For ja,p,i{s), we obtain the differential equation, 

3a,p,l{s) + {P + 3)ia,p+l,l(s) = QQ^^ ^ [p{(50 - 150^a,odd)s 

+ (180 - 780(5„,odd)s^ + (240 - 1650(5a,odd)s^ 
+ (150 - 18305«,odd)s^ + (45 - 1125(5«,odd)s' 
- 360(5c,,oddS - 45(5a,odd} 
+ {(46 - 150(5«,odd).s^ + (186 - 8105„,odd)s^ 
+ (330 - 1800(^«,odd)s^ + (360 - 2130(5«,odd)s^ 
+ (225 - 1395(5«,odd)s' - 450(5«,oddS - 45(5«,odd}] • (A-16) 
We assume that ja,p,i{s) is expressed as 

ja,p,lis) = P ■ Va,p{s) + '/'a,p(s), (A- 17) 

where i]a,p{s) and (t>a,p{s) satisfy eq. (A- 10). Under the above assumption, we substitute 
j'a,p,A^) ja,p+i,i{s) into the left-hand side of eq. (A-16) and obtain first-order linear dif- 
ferential equations for i]a,p{s) and <pa,p{s)- We solve the resulting equations and finally obtain 
ja,p,i{s)- Similarly, we assume that fa,p,2{s) is expressed as 

fa,p,2{s) = P^ ■ Ha,pis) + p ■ l^a,pis), (A-18) 

where /Xa,p(s) and i'a,p{s) satisfy eq. (A- 10). Repeating the procedure described above, we 
obtain 

fa,p,2is) = 3gQ^^ I ^y^, [{lip - 9 + (15p - 15)S^,om}s' 

+ {36p - 24 + (120p - 120)(5a,odd}s^ 

+ {90p - 60 + (240p - 240)5«,odd}s^ 

+ {120p - 150 + (180p - 180)(5„,odd}s^ 

+ {90p - 90 + {45p - 45)(5c,odd}s^] . (A- 19) 

To derive fa,p,3{s), we must add the differential equations for ya,p,o, Ua,p,o, ta,p,o, la,p,i, ka,p,i, 
ja,p,2 and fa,p,3{s) to eqs. (A-l)-(A-7). We present only the final result: 

^"-»<') = 45360('+.)--^e [«^"''' - + 

- (693p^ - 1953p + 1260)5„,odd}s^ 

+ {(756p2 - 1368p + 630) - (5544/ - 15624p + 10206)(5„,odd}s^ 
+ {(1890p2 - 2448^9 + 1008) - (19404^^ - 54684p + 36288)5a,odd}s^ 
+ {(2520/ - 2142p + 1260) - (36288/ - 106218p + 73458)(5a,odd}s^ 
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+ {(1890^2 + 3402p - 4914) - (39690j9^ - 120960p + 83916)5a,odd}s^ 
+ {(9450p - 24570) - (24570^^ - 75600p + 49140)(5a,odd}s'^ 
+ {(18900p - 24570) - (6615/ - 19845p + 11340)(5a,odd}s^ 
+ (5670^- 5670)8^1. (A-20) 
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